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Quantum dots pose a problem where one must confront three obstacles: randomness, interactions 
and finite size. Yet it is this confiuence that allows one to make some theoretical advances by 
invoking three theoretical tools: Random Matrix theory (RMT), the Renormalization Group (RG) 
and the l/N expansion. Here the reader is introduced to these techniques and shown how they 
may be combined to answer a set of questions pertaining to quantum dots. 
PACS 73.21, Tl.lOAy 
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I. INTRODUCTION 

This colloquium is based on a lecture entitled " Dots for 
Dummies" I have frequently given. The title was chosen, 
not to offend, but to keep experts in the audience from 
hijacking the lecture with minutiae while the intended 
goal was to introduce certain problems involving quan- 
tum dots to a broad audience not necessarily working in 
this subfield. The present article attempts to do the same 
for the general readership of this journal. To follow this 
colloquium you must be familiar with the rudiments of 
second quantization and Feynman diagrams. 

The emphasis of this article is idiosyncratic. I am par- 
tial to what I know best and like to talk about most: a 
frankly pedagogical tour of several useful techniques and 
their applications to this problem. There is not a compa- 
rable emphasis on phenomenology, for which I will direct 
you to other excellent sources. 

So let us begin by asking "What are some of the ques- 
tions in the world of quantum dots, what are some of the 
tools brought to bear in addressing them, and what are 
some of the answers?" 



For our purposes, the dot is an island of size L (in the 
nanoscale ) within which electrons are restricted to live. 

If the dot is very small and the level spacing large 
(compared to other scales like temperature and interac- 
tion strength) we may be interested in just a few quantum 
states of the dot. For example if the dot is essentially a 
two-state system, one may study it with the intention 
of using it as a qubit in quantum computation. Here 
one needs to understand in detail the level structure of 
a particular dot so we can manipulate (program) it in a 
controlled way. 

The focus here is on large dots with a very fine level 
spacing. These dots are closer to the bulk system of inter- 
acting fermions in that a large number of energy levels are 
in play. However the finite size and finite level spacing are 
essential features. These dots are very much like nuclei- 
finite systems of interacting fermions- with the main dif- 
ference that the confining potential is externally imposed 
and not internally generated. The hard-wall boundary of 
the dot is assumed to be sufficiently irregular so that clas- 
sical motion is chaotic at and around the Fermi energy. 
The dot is otherwise dirt-free and motion within is bal- 
listic. The system is assumed to be quantum-coherent 
across the length of the dot. 
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FIG. 1 A typical situation where leads bring in electrons 
which tunnel into and out of dot (opposite to the current 
shown by arrows) . A gate voltage Vg is the control parameter 
and V is the vanishingly small source-drain voltage difference. 
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FIG. 2 A caricature of conductance G versus gate voltage Vg. 
The peaks have varying heights, widths and spacings. 



The dot is juxtaposed with leads that act as source and 
drain for electrons. Electrons are allowed to tunnel in and 
out of the dot. A gate voltage Vg allows one to vary the 
dot energy relative to the Fermi energy of electrons in 
the leads as in Figure H]). A tiny voltage V is applied 
between the source and drain to drive the current, and 
the conductance G = I /V (in the limit of vanishing / 
and V) is measured as a function of the gate voltage Vg 
sketched in Fig ©.^ 

The theoretical challenge is to describ e the observed 



series of peak p o sitioii s and peak hei g hts (lAleiner et al. 
2OO2I : lAlhassi d'. '2000'. IChang et ali Il996l: iFolk et al. 



Guhr. MuUer-Groeling and Weidenmulleil . Il998 ; 



Jalabert. Stone and Alhassidl . 19921 : Merlinl . 



These are decided by the shape of the empty dot 
coupling to the leads, the number of electrons in it and 
the interactions between them. Now, no one disputes 
that given enough such information, every peak may be 
described in greatest detail since the underlying laws 
(nonrelativistic quantum mechanics and the Coulomb 
interaction) are thoroughly understood. However this 
is not the goal here. The goal is more akin to that 
set by Wigner in describing excited states of nuclei. A 
nucleus, like a dot, is a finite many-body system with 



strong interactions between its constituents. Wigner 
argued that while one could (in principle) describe the 
specific levels of a particular nucleus with its incredibly 
complicated Hamiltonian, it would more interesting to 
describe statistical properties of the spectrum. More 
specifically he suggested the following program. Consider 
a mathematical ensemble of Hamiltonians which have 
the same symmetries as the nuclear Hamiltonian (say, 
being real and hermitian) and have the same average 
level spacing. Calculate the statistics of level spacings 
in this ensemble. Compare this to the actual levels by 
averaging over members of the physical ensemble. ^ 

Note the difference with Statistical Me- 
chanics, as emphasized by Gu hr et al 
( Guhr. MuUer-G roeling and WeidenmuUerl . 



There we consider an ensemble of systems with the same 
Hamiltonian but different initial conditions, while here 
we consider an ensemble of Hamiltonians with the same 
symmetry properties as the given nuclear Hamiltonian. 

In the case of dots, the members of the physical en- 
semble contain dots with the same density of states. One 
can generate members of the ensemble by varying the 
number of electrons in one dot (which probes different 
energy regimes) or by manufacturing many similar dots 
of the same area. If dots of different sizes are used, one 
can rescale their energies so that the mean spacing is the 
same. The mathematical ensemble consists of Hamilto- 
nians of the same symmetry class (e.g., real) and same 
average level spacing. In cases where there is a magnetic 
field (and the ensemble contains complex unitary Hamil- 
tonians) one can also vary the field to move around the 
ensemble. (A minimum change in field may be needed to 
generate an independent member.) 

Here are two kinds of questions one could ask: 

• If I measure A, the difference in gate voltage be- 
tween successive peaks, what will be the probability 
distribution P(A)? 

• If I measure G, the maximum of each conduc- 
tance peak, what will be the probability distribu- 
tion P{G)1 

There are of course other issues that we could discuss, 
but will not, such as the width of each peak as a function 
of temperature or spin content. 

To get a feeling for the results, let us see how we would 
go about explaining the data, armed with just the theory 
of non-interacting electrons (whose spin will be initially 
ignored.) 

To this end let us consider a particular dot of some 
particular shape. We would naturally begin by solving 
the single-particle Schrodinger equation 



Ho(j)a{r) = £a(/>a(r) 



(1.1) 



There is a proportionality factor relating the actual gate voltage 
to the applied gate voltage which we have set equal to unity. 



^ It can sometimes happen that any large segment of the spectrum 
of almost every member of the physical ensemble exhibits these 
statistical trends, in which case we can look at just one member. 
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with (j)a{j^) = at the boundary. This would have to 
be done on a computer for a generic dot. Imagine that 
we have all the wave functions 4>a and energies £„. Let 
6i = £i+i — £i be the spacing between levels i and i + 1 
and let 5 stand for a generic one, as well as the average 
level spacing (which is the inverse density of states). 

In the non-interacting theory, it is easy to see that at 
a randomly chosen value of gate voltage Vg, there will 
be no conduction. This is because the Fermi energy of 
the electrons in the lead will typically lie between occu- 
pied and empty levels in the dot, so that transmission 
by the occupied levels is forbidden by the Pauli princi- 
ple and transmission by the empty levels is forbidden by 
energy conservation. However, when the Fermi energy of 
the leads equals an energy level of the dot there will be 
transmission. Note that during transmission the dot has 
the same energy with N or N-l-1 electrons. This feature 
survives even when interactions are included. 

In this free-particle model, the difference in the 

gate voltage Vg between peak i + 1 and peak i times — e, 
the charge of the electron, will equal i5i_i_|_i = Si+i — 
Ei, the difference in energy between the level that was 
responsible for peak i+1 and the one that was responsible 
for peak i. It follows that (suppressing the constant — e) 

P(A)-P(<5). (1.2) 

In other words, the distribution of spacings in Vg between 
successive peaks is the same as the distribution of level 
spacings in the dot. Consequently, to obtain the statisti- 
cal distribution P(A) theoretically, we need to solve for a 
large number of energy levels in one particular dot, record 
the spacings S, and repeat within the ensemble, which 
means here other dots with similar density of states but 
arbitrary shape. (In order to fold in the results from 
many dots and many energy ranges, it will be necessary 
to rescale energies so that S, the average level spacing, is 
fixed at some value.) 

As for the height of any conductance peak, it will be 
given by the product of the probabilities for hopping on 
to the dot at the left lead L, and hopping off on the right 
lead R. These are determined 4>a{L/R), the value near 
the leads of the wave function (f>a corresponding to the 
state which is responsible for the transmission. One could 
collect the statistics of these as well, from our numerical 
solutions to Eqn. (|l.ip . 

Now it turns out we can spare ourselves a lot of trou- 
ble in obtaining these probability distributions P{5) and 
P{G) by appealing to Random Matrix Theory or RMT, 
provided certain condition apply. So we begin with a 
crash course on RMT. 

The first step is to acquaint ourselves with Et = 
hvp/L^ the Thouless energy, where is the Fermi ve- 
locity. Evidently Et stands for the uncertainty in energy 
of an electron that traverses the dot ballistically at the 
Fermi velocity in a time r = L/vp- Thus if the dot is 
connected to big fat leads. 



single particle levels in this energy window will each con- 
tribute a maximum of /h to conductance G. Thus in 
this case of a dot connected strongly to leads, g will be 
the dimensional conductance, i.e., G measured in units 
of the quantum of conductance e^/ft.. . 

While this is one way to introduce Et, you could object 
that in the experiments we consider the leads are weakly 
coupled to the dot and the levels are sharp. You would 
be right, and it turns out we are interested in Et for the 
following different reason. 

If the dot has a generic shape with no conserved quan- 
tities except for energy, and the classical dynamics at the 
Fermi energy is chaotic, the statistics of energy levels ly- 
ing within a hand of width Et and of th e correspondin g 
wave functions, are determined by RMT jiMehtd . \l99i), 
given just general features like the average level spacing 
6 and symmetries of the Hamiltonian. This is like say- 
ing that the statistical properties of an isolated box of 
gas are determined by just gross features like volume, 
number of particles and energy. Such a statistical ap- 
proach to the energy levels of very complicated systems 
was first taken by Wigner in the case of nuclei, which 
are once again finite systems made up of a large number 
of particles subject to strong interactions that defy any 
direct treatment. Notice that to use RMT you just need 
to know the symmetry class (real or complex matrix el- 
ements etc.) but no further details including even the 
spatial dimensionality of the problem. 

The first RMT prediction is that P{5), the distribu- 
tion of exact energy level differences Si+i — Si (for states 
lying with a band of width Et), will be indistinguish- 
able from that of an ensemble of random matrices of the 
same symmetry class, which in our problem, consists of 
all real hermitian matrices with the same S, the average 
level spacing. (For a dot in a magnetic field, where time- 
reversal symmetry is broken, the ensemble will consist of 
complex hermitian matrices.) 

In other words, the results one person gets by painstak- 
ingly solving for the spectrum of an ensemble of dots and 
compiling the level statistics is same as that compiled by 
another person who picks real hermitian Hamiltonians 
out of a hat, diagonalizes them and plots their spacing 
distribution! Furthermore this universal distribution is 
known in analytic form as soon as we provide S. 

So, history repeats itself: when things get really bad, 
(the dynamics goes from integrable to chaotic) they be- 
come good again because statistical methods become ap- 
plicable. 

The second RMT result bears on wave functions. A 
common quantity of interest is the ensemble average (de- 
noted by (•••)) of products of wave functions such as 

mr)Mr')) (1.4) 

defined as follows. Pick one realization of the dot and 
number the states by energy, with a and /3 being two such 
numerical labels. Pick two points r and r' inside the dot. 
Find 0* (r)0^(r'). Now change the dot smoothly, staying 
within the ensemble. Since there is no level crossing in a 



4 




WOF basis (labeled by k) via the functions 
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FIG. 3 The Wheel-of- fortune states within a band of energy 
Et concentric with the Fermi circle. There are roughly g such 
states of mean momenta k centered on g equally spaced points 
on the Fermi circle. The WOF states are obtained by chop- 
ping off plane waves of the desired mean momentum at the 
edges of the dot. These states are very nearly orthonormal. 



chaotic dot, the labels a and P retain their integrity. Now 
recompute 0* (r)(/)^(r'). Find the average of such terms 
over the entire ensemble. This is what Eq. (II. 4|) means. 
Such correlat ions are of use in computing peak -height 
distributions ( Jalabert. Stone and AlhassidI Il99^ . 

In our problem we need the correlation of wave func- 
tions in momentum space rather than coordinate space 
and we shall use them to deal with electron-electron in- 
teractions rather than to calculate P{G). 

Let us begin by modifying the definition of "momen- 
tum space" as appropriate to the dot. 

Consider a circular Fermi sea in k space and a concen- 
tric annulus of width Et in energy. In the bulk, this re- 
gion contains an infinite number of k states. If we now go 
to the dot of size L, the best we can do vis-a-vis momen- 
tum is wave packets centered at some k and of width 1/ L 
in both directions. It is readily verified that we can form 
g such "Wheel-of-fortune " (WOF) states (as in Figure 
3l)) w ithin this annulus ( Murthv. Shankar and Mathud . 
2005] . Suppose we expand the g exact eigenstates (la- 
beled by a) of a dot within the Thouless band in this 



la) =5]0„(k)|k). 



t(k) as fol- 



(1.5) 



Then a typical RMT result invoked here is that as g oo 



(C(k)^Mk')) = 



'5kk"5a^ 

9 
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(1.6) 



where the < • • • > denote describe an average over an 
ensemble of similar dots. 

Note that Eqn. (|1.6p is the minimal correlation we 
must have: in each sample, a and k label two orthonor- 
mal bases, so that if we set k = k' and sum over k we 
must get Sa/s, sample by sample and hence on average as 
well. (The same goes for setting a — (3 and summing over 
them to get (5kk'-) Similar correlators exist for products 
of four wave functions and these have the form of Wick's 
theorem. 

Unlike the result on level spacings, which can be 
demonstrated analytically, this one on wave functions is 
an assumption we will make. A reader armed with the 
requisite determination and computing power is invited 
to confirm or contradict this assumption. All we can say 
is that several consequences of t his assumption have been 
verified in our numerical work ( Murthv et all l2004j ) . 

Before proceeding let me address a common question. 
Is there any reason to believe that the g exact eigen- 
states within Et can be expanded in terms o f the g WOF 
states? We (Murt hv. Shankar and Mathuil . [2009) have 
verified the following in a numerical study of a " billiard" , 
or dot. First we manufactured the g states of mean mo- 
mentum k by choosing g equally spaced points on the 
Fermi circle and then chopping off the plane waves of 
these momenta at the edges of the billiard. (We chose 
5 = 37 in our study.) Then we verified that these 
wavefunctions were orthonormal to an excellent accuracy. 
Next we asked how good a basis these functions formed 
for expanding the a states within Et- We found the ex- 
act eigenstate at the middle of the Thouless band, i.e, at 
the Fermi energy, retained more than 99.9% of its norm 
upon projection to the WOF basis. As we left the cen- 
ter of the band the overlap decreased and dropped to 
around 50% at the edges. Thus our results get more reli- 
able as we go deeper into the Thouless band centered at 
the Fermi energy. 

Let us turn to the data on P(A), the distribution of 
spacings in Vg between peaks. In the noninteracting case 
we saw that P(A) = P{S), and then we saw that P{6) 
was given by RMT. 

So let us test the RMT results for level spacings P{S), 
against -P(A), the measured distribution of spacings in 
Vg between successive peaks. These will match only if 
the assumptions of non- interacting electrons is correct. 

The test fails. It is found that the peak spacings are 
an order of magnitude larger than level spacings in the 
dot. The reason is of course electron-electron interac- 
tions. Even if we ignore the full quantum mechanical 
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treatment of interactions, we need to acknowledge the 
classical fact that when we add an electron to a dot, it 
experiences Coulomb repulsion for the other occupants. 
Thus the energy cost of adding an electron is not just 
that of going to the next empty level, it is the additional 
repulsive energy. This capacitivc charging energy is pro- 
portional to and we must add a term uqN^ to the 
second quantized Hamiltonian. By studying the data one 
can come up with a good fit to uq. 

However we can anticipate that this alone will not work 
in the presence of spin, even in the non-interacting case. 
With spin present, each single particle orbital can be dou- 
bly occupied. Thus adding an electron will (will not) cost 
us a single particle energy level difference if it brings the 
total occupancy to an odd (even) integer. Combined with 
the above mentioned charging energy, we expect the peak 
spacing to have a bi-modal distribution as N varies from 
odd to even. 

This is however not seen. The reason is based on the 
exchange interaction. While it makes sense to put two 
electrons in each orbital if one is minimizing kinetic en- 
ergy, this may be a bad idea in terms of potential energy: 
the electrons in any given orbital would have opposite 
spins and hence the Pauli principle would not keep them 
away from each other, thereby raising the Coulomb repul- 
sion. To minimize the potential energy we must occupy 
each orbital once and place the electrons in the same spin 
state so that the Pauli principle can keep them apart in 
space. To let the system decide what is best we must give 
it an incentive for higher spins by adding a term —JS^, 
where S is the total spin. 

Thus we are led to the following 

S second-quant ized) " Universal" Hamiltonian: 

Aleiner et al\ . 
1998 
2000 



i2ooa 



Andreev and Kamenev 



Baranger. Ullmo and Glazman 



Browner. Oreg and Halperinl. 



1998; 



Kurland.Aleiner and Altshuleil . l2000l: lOrege^ aLl . l200lf ) 



H, 



(1.7) 



where ipa destroys an electron in a state a. A third 
possible interaction term pertaining to superconducting 
fluctuations has been dropped on the grounds that it is 
unimportant for the dots in question. 

Some proponents of the universal Hamiltonian give the 
following argument for why no other interactions need be 
considered. Suppose we take any other familiar interac- 
tion and transcribe it to the exact basis. Sums of prod- 
ucts of random wavefunctions 4>a will appear and lead to 
terms with wildly fluctuating phases. These terms can be 
shown to have zero ensemble average as g — > oo. Since 
deviations from the zero average will be down by 1 /g we 
can drop them at large g. By contrast the two terms 
kept (which commute with Hq) are free of oscillations 
and survive ensemble averaging. 

While the success of Hjj in explaining a lot of data 
is unquestioned, the accompanying arguments are not 
persuasive. In particular, ensemble averages should be 



performed not on the Hamiltonian but on calculated ob- 
servables. It is also not clear that a group of terms should 
be dropped because they are small, since they could ul- 
timately prove important to physics at very low energies 
(low, even within the already tiny band of width Et)- 

Now there is a tried and tested way for determining the 
relative importance of possible interactions in the low en- 
ergy limit. It is the Rcnormalization Group (RG). Not 
only can it tell us how important any given interactions 
is in the low energy sector, if there are competing inter- 
actions, it can provide an unbiased answer in w hich they 
are a ll allowed to fight it out. (As explained in ( Shankan . 
Il994[ ). The Luttinger Liquid in d = 1 is an example in 
which Superconductivity and Charge Density Wave for- 
mation compete and actually annul each other while in 
d — 2 the latter wins at half-filling on a square lattice.) 

It is therefore natural to ask if the RG applicable in 
this problem and if so what its verdict is. 



II. THE RENORMALIZATION GROUP - THE TEN CENT 
TOUR 

In order not to leave behind readers unfamiliar with the 
RG the following lightening review is provided. Experts 
can skip to the next section. 



A. What is the RG? 

Imagine that you have some problem in the form of a 
partition function 



Z{a, b) — j dx J dye 

dx ( dye'^^^'^y'''''''-^ 



(2.1) 
(2.2) 



where a and b (along with other such possible terms) are 
constant parameters and S is called the action. (In clas- 
sical statistical mechanics it would be called the energy.) 

The parameters b, c etc., are called couplings and the 
monomials they multiply are called interactions. The 
x^ term is called the kinetic or free-field term, and o, the 
coefficient of the kinetic term, has no name and is usually 
set equal to ^ by rescaling x. 

The average of any f(x,y) is given by 



/da;/ dy/(a;, y)e~''(^^+f^^e~^(^+2^'* 

JdxJdye-^(-'+y')e-b(^+v^^ ' ^^'^^ 

Suppose we are interested only in functions of just x. 
In this case 



6 



/rfx/(x)/dye-°("'+^')-''("+^)' 
f dx/(a;)e-"'^'-^'^''+- 



(/(^)) 



where 



-a' x'^ —b' x'^ .. 



dye 



(2.4) 



(2.5) 



defines the parameters a', 6'. ..etc. In the general case we 
would like to define an effective action S'(x; a', 6', c'..) by 

g-S'(x;a',6'.c'..) ^ j j^y^-S{x,v;a,b,c,..) ^^.Q) 

and the corresponding partition function 



Z{a',b'...) = / dxe-^'^"-'^''''^' -) 



(2.7) 



which define the effective theory for x. These parameters 
a', h' etc., will reproduce exactly the same averages for x 
as the original ones a, 5, c. did in the presence of y. This 
evolution of parameters with the elimination of uninter- 
esting degrees of freedom, is called renormalization. It 
has nothing to do with infinities; you just saw it happen 
in a problem with just two variables. 

Note that even though y does not appear in the effec- 
tive theory, its effect has been fully incorporated in the 
process of integrating it out to generate the renormalized 
parameters. We do not say " We are not interested in 
y, so we will set it equal to zero everywhere it appears" , 
instead we said, "What theory, involving just a;, will give 
the same answers as the original theory that involved x 
and y?" 

The term "Group" arises as follows. Let us say y is a 
short hand for many variables as is always the case in real 
life. If we eliminate one of them, say yi which leads to 
the renormalization (a, 6, c, ..) — » (a', 6', c'..) and then we 
eliminate y2 so that now (a', 6', c', .■) (a",6",c"..) the 
net result is equivalent to a single renormalization process 
in which {a,b,c, ..) (a",6",c"..) under the elimination 
of 2/1 and y2- (The RG is not a real group since we cannot 
define a unique inverse.) 



B. How is mode elimination actually done? 

Notice that to get the effective theory we need to do a 
non-gaussian integral. This can only be done perturba- 
tively. At the simplest Tree Level, we simply drop y and 
find b' = b. 

In other words, at tree level, the effective action is 
found by simply setting the unwanted variables to zero 
in the original action. 



At higher orders, we bring down the non-quadratic ex- 
ponential and integrate in y term by term and generate 
effective interactions for x. 

Here is how it is done in our illustrative example. 

2 4 r^]„,p-ay^ -b{ixy^+ix''y+6x^y^+y^} 



Zoia) 



e Zo{a){e 



-ax--bx'^r^ I „\ I „-b(Axy^+Ax^y+6x^y'^+y'^)\^ ^ (2 8) 



where we have multiplied and divided by ZQ[a), 

Zo(a) = [ dye-^y^ (2.9) 



the partition function for a Gaussian action e and 
where 



-b(Axy^+ix^y+&x^y'^+y'^) 



(2.10) 



stands for the average of the exponential with respect 
to the partition function Zo{a). Since ZQ{a) is indepen- 
dent of a;, we will simply ignore it in the effective action 
S'{x; a' ,b' ..), without altering any absolute probability. 
As for Eq. (|2.10|) we invoke the result valid for averages 
over Gaussian actions: 



(2.11) 



This is called the cumulant expansion and we see that in 
our example where V = b{Axy^ + Ax^y + Qx'^y^ -\- y*), the 
exponent is a power series in 6 for contributions to S". 
To the leading order in 6 we have 

S'{x-a\b'..) = ax^ + bx'^ + {biAxy^^ -^Ax^y + Qx^y^ + y^)) 
= ax^ + bx^ + 6bx^(y'^) + (y^) 

= ax'' + bx^ + 6b— + ^ (2.12) 

because (y) = {y^) = 0, {y') = ^ and (y'*) = Thus 
we have our first RG result for renormalization: 



3b 
a 



(2.13) 



to leading order in b. Note that we do not care about 
^-independent constants like For reader who know 
Feynman diagrams, the power series in b can be identified 
with the Feynman graphical expansion for a (j)'^ theory 
(since we kept up to quartic terms in the action.) In 
the case of the actual field theory itself, when the 
momentum cut-off is reduced from A to A/s (where s > 
1), X will denote low- momentum modes < fc < A/s and 
y the high- momentum modes A/s < fc < A, and the loop 
integration will be over the range A/s < fc < A. 
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C. Why do the RG? 

Why do we do this? Because the ultimate effect of any 
couphng on the fate of x (the variable we care about) is 
not so apparent when y (the variable in which we have no 
direct interest) is around, but surfaces to the top only as 
we zero in on x. For example, we are going to consider a 
problem in which x stands for many low-energy variables 
and y for many high energy variables. As we integrate 
out high energy variables and zoom in on the low energy 
sector, an initially tiny coupling can grow in size, or an 
initially impressive one diminish into oblivion. 

This notion can be made more precise as follows. Con- 
sider the gaussian model in which we have just a ^ Q. 
We have seen that this value does not change as y is 
eliminated since x and y do not talk to each other. This 
is an example of a fixed point of the RG since the cou- 
pling that goes in comes out unchanged by elimination of 
unwanted variables. ^ Now turn on new couplings or " in- 
teractions" (corresponding to higher powers of x, y etc.) 
with coefficients 6, c and so on. Let a', h' etc., be the 
new couplings after y is eliminated. Let us in addition 
rescale x so that x"^ has the same coefficient as before 
i.e., a' = a = ^say. ^ Any of the renormalized couplings 
(still called b',c',..), which are are bigger than the ini- 
tial ones, 6, c, d.., are called relevant while those that are 
smaller are called irrelevant. This is because in reality 
y stands for many variables, and as they are eliminated 
one by one, the relevant coefficients will keep growing 
with each elimination and the irrelevant ones will keep 
shrinking and ultimately disappear. If a coupling neither 
grows not shrinks it is called marginal. Thus the RG will 
tell us which couplings really matter in the low energy 
limit which controls the nature of the ground state and 
its low lying excitations. 

There is another excellent reason for using the RG, and 
that is to understand the phenomenon of universality in 
critical phenomena, an area I will not enter. 

Later we will apply these methods to quantum dots. 
For the reader who may understandably get restless in 
the interim, here is a glimpse how we will use RG in the 
case of dots. First we will begin with all single-particle 
states in the Hilbert space of dot. Then we will zero 
in on states lying within a narrow band of states within 
an energy El the Fermi energy Ep- (AH one needs is 
that El << Ep.) In other words, at this point the vari- 
ables we called x {y ) are states lying inside (outside) 



this band. For the clean system in t he bulk, the res ult 



^ There can be more complicated fixed points where this happens 
despite the fact that the x's and j/'s talk to each other, and the 
integrals do not factorize. While we will not encounter them in 
our problem, the proposed strategy applies there as well. 

* We like to keep a fixed because what really matters is the relative 
size of a and the interaction terms. For example if it turns out 
that V > b but also that a' > a, it is not clear the renormalized 
theory has stronger interactions. Keeping a fixed allows us to 
compare apples to apples. 



of such a mode elimination is known ( ShankaJ . Il994^ to 
lead to the result that the system is described by an in- 
finite number of Landau Fermi Liquid parameters Urm 
where m is an integer. This process will be discussed 
in abridged form. The Universal hamiltonian contains 
just the uo term. ^ We would like RG to tell us if and 
when we can ignore all other Wm's. To this end we will 
begin with a theory where only states within Ep of Ep 
are kept and the starting hamiltonian has all Landau in- 
teractions written in terms of the dot cigcnfunctions <f>a ■ 
Then we will ask what happens to them as we eliminate 
states within Et, getting even closer to Ep. To see what 
happens when we do this (and how we do this) you need 
to read on! 



III. THE PROBLEM OF INTERACTING FERMIONS 

We will now apply RG to a system of nonrelativistic 
spinless fermions of mass m and momentum K in two 
space dimensions. The single-particle Hamiltonian is 



H 



2m 



(3.1) 



where the chemical potential fi is introduced to make sure 
we have a finite density of particles in the ground state: 
all levels within the Fermi surface, a circle defined by 



2m 



(3.2) 



are now occupied since occupying these levels lowers the 
ground-state energy. 

In second-quantization we are thus starting with the 
Hamiltonian 

H, = j dPKij\K) (^^^^^) V^K) (3.3) 

where "(/"^(K) creates a fermion of momentum K. No- 
tice that this system has gapless excitations above the 
ground state. You can take an electron just below the 
Fermi surface and move it just above, and this costs as 
little energy as you please. Such a system will carry a dc 
current in response to a dc voltage. An important ques- 
tion one asks is if this will be true when interactions are 
turned on. For example the system could develop a gap 
and become an insulator. Or it could become a super- 
conductor. How do we decide what the fermionic ground 
state and low energy excitations will be, short of solving 
the problem? 

In the noninteracting limit the ground state is simple: 
it is a single Slater determinant (antisymmetrized wave 



There are actually two sets of u's, one for charge and one for 
spin. The two no's are the charging and exchange interactions 
UQ and Jq. 
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FIG. 4 The low energy region for nonrelativistic fermions 
lies within the annulus concentric with the Fermi circle. It 
extends A in momentum and -El in energy from the Fermi 
circle. In units where the Fermi velocity vf is chosen to be 
unity, the two are equal. 



function) with all momentum states below Kp occupied. 
If we turn on say, quartic interactions, the new ground 
state will have pieces in which two fermions from below 
Kp have been moved to two fermions above Kp keeping 
the total momentum at zero. Each such piece will come 
with a numerator proportional to the interaction strength 
and an energy denominator equal to the cost of creating 
this particle-hole pair. Clearly fermions deep in the sea 
will not take part in this kind of process with any serious 
likelihood as long as the interactions to be added are 
weak. Likewise if states far above Kp were removed no 
one would be the wiser. Thus the fate of the system 
is going to be decided by states near the Fermi energy. 
This is the great difference between this problem and 
the usual ones in relativistic field theory and statistical 
mechanics. Whereas in the latter examples low energy 
means small momentum, here it means small deviations 
from the Fermi surface. Whereas in these older problems 
we zero in on the origin in momentum space, here we 
zero in on a surface. The low energy region is shown in 
Figure IH 

We are now going to add interactions and see what 
they can do. This is going to involve further reduction 
of th e cut-off. Here we join the discussion of ( Shankail . 
I1994D . 

Let us begin then with a momentum band of width A 
on either side of the Fermi surface. (In terms of energy, 
the band has a width E^, where L stands for Landau, 
for reasons that will follow. ) 

Let us first learn how to do RG for noninteracting 
fermions. To apply our methods we need to cast the 
problem in the form of a pa th integral. Following any 
number of sources, say (| Shankail . Il994l ) if you want to 
one-stop shopping, we obtain the following expression for 
the partition function of free fermions: 



,So 



(3.4) 



where 
and 



2m 



[dijj(h/J\ = J|dV'(w,K)dV'(w,K) 



(3.5) 



(3.6) 



K 



and and ijj are called Grassmann variables. They are 
really weird objects one gets to love after some familiar- 
ity. There are some rules for doing integrals over them. 
For now, I suggest you note only that (i) ^/;(a;,K) and 
ijj{uj, K) are defined for each to and k in the annulus (ii) 
Z is a product of Gaussian integrals, with the variables 
at each (u;,K) not coupled with those at anot her. The 
dedic ated reader can learn more from Ref. ( Shankail . 
I1994D . 

We now adapt this general expression to a thin annulus 
to obtain 



Zo = j [d^d^ 



(3.7) 



where 
So 



27r poc />A 

dd duj I dkip{iuj — vp k)'4i. (3. 
J-oo J-A 



We have approximated as follows: 
K'-Kl ^Kp 

— hj — up n- 



2m 



m 



(3.9) 



where k — K — Kp and vp is the Fermi velocity, hereafter 
set equal to unity. Thus A can be viewed as a momentum 
or energy cut-off measured from the Fermi circle. We 
have also replaced KdK by Kpdk and absorbed Kp in 
Tp and tp. It will be seen that neglecting k in relation to 
Kp is irrelevant in the technical sense. 

Let us now perform mode elimination and reduce the 
cut-off by a factor s. Since this is a gaussian integral, 
mode elimination just leads to a multiplicative constant 
we are not interested in. So the result is just the same ac- 
tion as above, but with |fc| < A/s. Consider the following 
additional transformations: 



iu;',k') - s{u;,k) 



{tp'{uj',k'),'ilj {uj',k')) 



uj' k' 



(3.10) 



(3.11) 



When you perform the highly recommended exercise 
of making this change of variables, you will find that the 
action and the phase space all return to their old values. 
So what? Recall that our plan is to evaluate the role 
of quartic interactions in low energy physics as we do 
mode elimination. Now what really matters is not the 
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absolute size of the quartic term, but its size relative to 
the quadratic term. Keeping the quadratic term identical 
before and after the RG action makes the comparison 
easy: if the quartic coupling grows, it is relevant; if it 
decreases, it is irrelevant, and if it stays the same it is 
marginal. 

Let us now turn on a generic four-Fermi interaction in 
path-integral form: 

Si^ J ^(4)V^(3)V(2)V'(l)^i(4, 3, 2, 1) (3.12) 

where / is a shorthand: 

^Yl d0, dh / dto, (3.13) 

At the tree level, we simply keep the modes within 
the new cut-off, rescale fields, frequencies and momenta 
, and read off the new coupling, a highly recommended 
exercise. We find 




FIG. 5 Kinematical reason why momenta are individually 
conserved up to a permutation. 



origin, to which the low energy region collapses. Here the 
low energy manifold is a 2-sphere no matter how small A 
is and uq will inevitably have dependence on the angles 
on the Fermi surface: 



u'{k\uj',e) = 




(3.14) 



This is the evolution of the coupling Junction. To deal 
with coupling constants with which we are more familiar, 
we expand the functions in a Taylor series (schematic) 

u = uq + kui + fc^U2--- (3.15) 

where k stands for all the fc's and w's and of course the 
9-s are fixed. An expansion of this kind is possible since 
couplings in the action are nonsingular in a problem with 
short range interactions. If we now make such an expan- 
sion and compare coefficients in Eqn. p.l4p . we find 
that uq is marginal and the rest are irrelevant, as is any 
coupling of more than four fields. For example 



Therefore in this theory we are going to get coupling 
functions and not a few coupling constants. 

Let us analyze this function. Momentum conservation 
should allow us to eliminate one angle. Actually it allows 
us more because of the fact that these momenta do not 
come form the entire plane, but a very thin annulus near 
Kp- Look at Figure m Assuming that the cutoff has 
been reduced to the thickness of the circle in the figure, 
it is clear that if two points Ki and K2 are chosen from 
it to represent the incoming lines in a quartic coupling, 
the outgoing ones are forced to be equal to them (not in 
their sum, but individually) up to a permutation, which 
is irrelevant for spinless fermions. Thus wc have in the 
end just one function of two angles, and by rotational 
invariance, their difference: 



Ui 



^ (3.16) 



Suppose we start with some cut-off Aq and an initial 
coupling, we shall call Wn(0). Let us now reduce A con- 
tinuously as per A = Aqc^* = A/s so that the change 
in Un is continuous in t. The relation u,j = u„/s", Eq. 
p.l6p . can be written as u„(i) = u„(0)e~"*. 

Sometimes this result is rewritten in terms of the /?- 
function; 



dUn{t) 

P[t) = — 7— = -nun(t) 
at 



where t = Ins. (3.17) 



If you consider the scalar field theory in four dimen- 
sions you find again an equation like Ea. (|3.15p with the 
difference that k is measured from the origin and not the 
Fermi surface. Thus all we have left to worry about is one 
number mq, the first term in the Taylor series about the 



%)^F{e^-e2) = u{e). 



(3.18) 



About for ty years ago L andau came to the very same 
conclusion! Landaul Il956t ). that a Fermi system at low 
energies would be described by one function defined on 
the Fermi surface. He did this without the benefit of the 
RG and for that reason, some of the leaps were hard to 
understand. 

Since u is marginal at tree level, (gets neither big nor 
small as the RG process is iterated and modes are elim- 
inated) we have to go to higher orders in perturbation 
theory to break the tie, to see if it ends up being relevant 
or irrelevant. The answer, given without proof, is that 
it remains marginal to all orders ( Shankarl Il994l ). This 



is for strange kinematical reasons. A brief review of how 
the higher order calculation will be done will be given 
when we come to dots. 
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Often one writes 



u{0) = Urn COs{md) 



(3.19) 



where the Landau parameters. The corresponding 

interaction^ is, in schematic form with radial integrals 
over k suppressed, 



Hl=J2 n{0)n{e')u„, cos {m{e - 0')) 



(3.20) 



where n{9) is the number density at angle 9 on the Fermi 
circle. 

Note that if uq is the only non- vanishing Landau cou- 
pling, we get the uqN'^ interaction of Hij. If we include 
spin, there are spin density-spin density interactions and 
Jo corresponds to keeping just the zeroth harmonic. So 
Hjj amounts to keeping just the lowest harmonic in the 
Landau expansion. 

We need to ask when and why m > terms can be 
ignored. 



IV. RG MEETS DOTS 

The first crucial st ep towards this goal was t aken by 
Murthy and Mathur (jMurthv and Mathuil . r2002f ) . Their 
strategy was as follows. 



• Step 1: Use the clean system RG described above 
to learn that at low energies the important interac- 
tions are the Landau interactions (of which uq and 
Jo are a subset). 

• Step 2: Start at the energy scale Et and switch 
to the exact basis states of the chaotic dot, writing 
the kinetic term and all the Landau terms in this 
basis. Run the RG by eliminating exact energy 
eigenstates within Et- See what happens to all 
Landau terms, do they grow, fall or stay the same? 

I will discuss some of the subtleties here, referring you 
to (jMurthv et all l2004 ) for more details. 

A common concern is to ask if we can ignore the walls 
of the dot and use momentum states at high energies. 
After all, momentum is not conserved in the dot and what 
meaning is there to using the Landau interaction written 
in the momentum basis? The point is this. Gonsider the 
collision of two particles in a box. As long as the collision 
takes place quickly (in terms of the time to cross the box) 
the incoming and outgoing particles can be labeled by 



momenta and this label will be useful (conserved) during 
the collision. ^ 

However this only brings us down to Ei^ , the energy be- 
low which Landau theory works. Although << Ep^ it 
is still finite for an infinite system while Et — 1/ Lis even 
smaller for large L. In this no man's land between El 
and Et^ the flow of couplings is intractable and landau 
parameters could have evolved from their bulk values. So 
what we are saying is this: the universal Hamiltonian has 
just two of the infinite Landau parameters (moj ^^o) in it- 
Let us put in the rest of them with any size, and see what 
their fate is under further mode elimination. 

While this looks like a reasonable plan, it is not clear 
how it is going to be executed. There are at least two 
obvious problems. To understand them you must un- 
derstand in more detail how the RG is done at the one 
loop level. To illustrate this we use the more familiar 
scalar theory in d = 4 with the UQ(f>^ interaction, where 
uo, hereafter called just u, is marginal and momentum- 
independent. Suppose we want the scattering amplitude 
of four scalar particles. The answer will be given by a 
perturbation series, depicted in Fig. ([6]). The second 
term is an integral over the loop momentum K that 
goes from < < A. The integrand is the prod- 
uct of the two propagators, which go as each at 
large K. (We have assumed the external momenta are 
all zero and neglected two other loop diagrams that be- 
have the same way.) The left hand side corresponds to 
a physical process and cannot depend on the cut-off A. 
It follows u must acquire A-dependence so as to make 
the right hand side A- independent. We can find u{K) as 
follows. Suppose A is reduced by \dA\. The loop will 
change by an amount equal to the integral in the region 
A — \dK\ < K < A, which is now missing. This miss- 
ing part must then come from the first term to keep the 
physical scattering amplitude fixed. A simple calculation 
shows that up to numerical factors 



so that 



, 2MAI 
A 



(4.1) 



(4.2) 



In summary, the change in u(A) is given by the loop 
terms with internal lines restricted to the modes being 
eliminated. 

We run into the following problems if we try to do the 
mode elimination for the dot. 

First, knowledge of the exact eigenfunctions is needed 
to even write down the Landau interaction (the analog 



While this form of will suffice for the dot, in the clean bulk 
system one must allow for small non-forward scattering since A 
is finite. Despite their small measure these terms are crucial 
because the Fermi liquid has a very singular response at the 
Fermi surface. 



There will be a parametrically small number of cases ( ratio of 
interaction range to dot size) where collisions take place near the 
walls, when this will be wrong. 
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>4- 



FIG. 6 Schematic of scattering amplitude in i^* theory. The 
sohd dot is the full answer and the empty one is the coupling 
constant that goes in. The loop is the first of many that enter. 



of u in the 0'* example) in the disordered basis: 



kk' 



X(0^(k')05(k)-0^(k)05(k')) 



(4.3) 



where k and k' take g possible values, 9 and 0' are the 
angles associated with k and k' and the interaction has 
been antisymmetrized to mate with the four-Fermi oper- 
ators ip'l'tp'l^'tp^'tps- 

Next, additional information is needed on energy lev- 
els to do the higher order (loop) calculation since the 
propagator for state a is (iu — Eq)"^. Remarkably it is 
possible to overcome ignorance of specific energy levels 
and wavefunctions. Here are some, but not all details of 
how this comes about. 

First consider a specific realization. We want to re- 
duce the cut-off by summing over some high energy-states 
(y's). If we write down expression for the one loop flow, 
four-fold products of the unknown wave functions appear 
at each vertex. At the left vertex, two of these wave func- 
tions correspond to external lines, which are fixed (and 
below the new cut-off, i.e., these are the x variables) while 
the other two correspond states to be eliminated and thus 
summed over. A similar thing happens at the right ver- 
tex. In addition there are the propagators for the two 
lines dependent on the single-particle energies. 

Thus what we have here is a product of four wave- 
functions and an energy denominator summed over states 
being eliminated. This sum will of course vary from dot 
to dot. Suppose we ignore this variation and replace 
the sum by its ensemble average. (The product of wave- 
functions and the energy denominators can be averaged 
independently since to leading order in 1/g, RMT does 
not couple them.) The entire average will then be just 
a function of just g as in Eq. (|1.4p . But what about 
sample -to-sample variations? Her e we invoke the nice 
result (|Murthv and Mathurj . |2002[) that that deviations 
from the average are down by an extra power of g. Thus 
the flow is self- averaging, and every dot will have the 
same flow as q ^ (x. ^ 



FIG. 7 The flow of the coupling Um for any m 0. The 
origin of this flow is the universal Hamiltonian where every 
Um = except uq which can have any value. All points to the 
right of u* flow to this. We ask what happens if we begin to 
the left. 



Using such a device, these authors found the remark- 
able result that the renormalized Va^^s is itself equivalent 
to a Landau interaction but with renormalized values of 
Um flowing as per 



dUri 



dt 



cu„ 



TO 7^ 



(4.4) 



where c is independent of m and of order unity. 

Note that uq does not flow and that j ust as in the BCS 
flow of the clean system ( ShankaJ . ll994l ). different to's do 
not mix to this order. If spin were included Jq wouldn't 
flow either. This lack of flow is due to the fact these co- 
efhcients multiply operators that commute with Hq, the 
non-interacting part of the Hamiltonian, and therefore 
have no quantum fluctuations. 

The flow (for to 7^ 0) implies that all positive Um's flow 
to zero, as do negative ones with Um > u* , the fixed point 
of the flow. Thus all points to the right of it* flow to Hu , 
as shown in Fig. [T] If a modest amount of Um, to > of 
either sign is turned on at the beginning (when A = Et): 
it will renormalize to zero as we go down in energy. On 
the other hand if we begin to the left of u* , we run off to 
large negative values. 

The universal Hamiltonian is thus an RG fixed point 
with a domain of attraction of order unity. To me this 
is the most satisfactory explanation of its success. 



V. THE l/N APPROXIMATION 

So far we have used the RG to understand the low en- 
ergy behavior of the quantum dot. We found that if we 
begin our analysis with states within Et of the Fermi 



Experts should note that we are not averaging the /3-function, 
it is self-averaging. There may also be concern that there will 



not be enough states within dK to allow the averaging to work. 
This too can be handled: there is a way to define the /3-function 
dShankarl . I1994I) in which the loop is first summed over all states 
inside A and then the A derivative is taken. A simpler illustration 
of self-averaging will be given later when we extend this result 
to all orders. 
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energy Ep and all possible Landau interactions Um, and 
slowly reduce the energy cut-ofF, we end up with the uni- 
versal Hamiltonian as the low energy fixed point for all 
positive Urn's and negative u^s that are are not more 
negative than some critical value From this knowl- 
edge at very low energies we can conclude in particular 
that the ground state (lowest of all energies!) is deter- 
mined by Hu for Um > and that at Um = the sys- 
tem undergoes a second order phase transition, for that 
is what a zero of the beta function or fixed point signifies. 
The RG does not tell us much about the nature of the 
new phase, other than that it will be dominated by the 
Um that runs off to strong coupling. 

While this is nice, there is clearly room for improve- 
ment. In particular, if the theory could be solved exactly, 
we would know even more. We would know if the phase 
transition at u,„ of order unity, found in the perturbative, 
one-loop calculation, really does occur. We would also 
know to what state the system is driven after the transi- 
tion. In other words, if we can solve a problem exactly, 
we do not need the RG. We could, after obtaining the 
solution with some cut-off, easily ask how the coupling is 
to be modified if the cut-off is then reduced. We would 
do this by computing a physical quantity P(u(A), A) as 
a function of the input cut-off A and input coupling u(A) 
and then set dP/dA = to find du/dA. This would be 
the exact beta function that one could compare to any 
perturbative result. 

It turns out we can do all this in our problem of dots. 
Now, any controlled calculation is based on some small 
quantity. If the small quantity vanishes the controlled 
calculation gives exact results. The most common exam- 
ple is the coupling constant itself, and the exact solution 
that emerges in the limit of zero coupling is a free the- 
ory. However, there are other cases exhibiting nontrivial 
behavior, where the coupling is not small but something 
else is. For our dots the small parameter turns to be 1/g 
and the exact solution emerges in the lim it 5 ^ oo as 
Murthy and I (jMurthv and Shankail . l2003l ) found. Note 
that the limit g — > cxd corresponds to the limit of infinite 
dot size. (Unfortunately this does not mean we under- 
stand bulk physics, since in this large dot the Thouless 
band, over which we have control using RMT, shrinks to 
zero as i/L.) What we assume is that what happens in 
this infinite dot will seen also in large but finite dots, just 
as one assumes that genuine phase transitions, which are 
allowed only in infinite systems will be well mimicked in 
large but finite systems. We could show that in this limit 
there is indeed a transition at a negative coupling of or- 
der unity and that to the left of it is a symmetry-broken 
phase which can be analyzed in some detail. 

We showed that one did not have to rely on RG or 
perturbation theory in powers of u. Instead the theory 
could be solved by saddle point methods for any u thanks 
to the smallness of l/g. The trick was to use a variant of 
the so called 1 /N expansion. So let us first get acquainted 
with the 1/N ex pansion. I give only a few details here, 
referring you to |Murthv a/.l . l2004) . 



Let us begin with the common situation wherein we 
expand the answer in a power series in the coupling u, 
via Feynman diagrams of increasing complexity. This 
method works provided the answer can be expanded in 
a series at u = and we limit ourselves to weak cou- 
pling (since we can typically compute only to some small 
order in u). There are however cases where we need to 
go to all orders in u or pick up essential singularities be- 
fore the right physics can be found. Amazingly some of 
these problems are tractable thanks to the 1/N expan- 
sion. Here is a brief survey of that trick. 

Suppose we have a theory of fermions with an internal 
isospin or flavor label that runs over N values. (The 
method works for bosons as well.) Let the theory be 
defined by the following schematic path integral: 



(5.1) 



where D stands for the quadratic kinetic energy term, 
and the sum over flavors or integral over space-time is 
suppressed, so that for example 



N 



dvdO Tpiipi 



The factor of 1/A^ in the quartic term is to ensure that 
the single sum over the flavor index in the kinetic term 
has a chance against the double sum in the quartic term. 
In a theory with Dirac fermions D = <^ while in a non- 
relativistic problem it could be d / dr — Ea- 

If we now introduce a Hubbard-Stratonovic field a we 
can rewrite Z as 



(5.2) 



the correctness of which is readily verified by doing the 
Gaussian integral over a. If we now use the fact that for 
each flavor 



= dei(i? + a) =e^'-'"(^+") (5.3) 

and that each flavor gives the same determinant, we ob- 
tain 

Z = j [d^d^ dae^{o+a)i,-No-l2u (5 4) 



It is now clear that in the limit N 00, we can do 
the fj integral by saddle point. At large and finite N, we 
can compute corrections in a series in 1/7V. However, the 
complete and exact dependence on u is obtained to each 
order in 1 /TV. Often just the leading term at = 00 cap- 
tures all the novel physics. A celebra ted example is the 
1 + 1-dimensional Gross-Neveu model ( Gross and Neveul . 
I1974D . where there is a coupling constant u and an in- 
ternal 0{N) isovector index that runs from 1 to N. In 
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the limit N ^ oo the saddle point method allows one to 
show that the a field spontaneously develops a non-zero 
average (which translates into a mass for the fermion) 
that goes as e~^/", a result clearly outside the reach of 
traditional perturbation theory. Corrections around the 
saddle point in powers of 1/A^ do not change the main 
features quoted above. In summary, the small parameter 
that makes a calculation possible is not u but 1/A^. 

We are going to do the same thing here, with 1/g being 
the small parameter. You may object that since in our 
problem D = d/dr — the g different fermions are 
not related by symmetry, and the appearance of a large 
number (the analog of N) in front of the Tr In is by no 
means assured. However, we found that if we went ahead 
and evaluated the Tr In order-by-order in a and exploited 
self-averaging as in the one-loop flow, a large number (g^) 
does indeed appear in front of the action, playing the role 
of A''. Here is a glimpse of how this happens. 

Let us first consider just one Landau term with cou- 
pling Umi which will simply be called u: 



Z 



[dixliP] dae'^(^+<-)V'-'-V2K, (5 6) 



Next let us make the change tr — > gti in Eqn. (|5.6p . 
While this trivially brings a in front of the cr^/2u term, 
what is nontrivial is that every term in the Tr In ( devel- 
oped in powers of a) will also go as g^ to leading order. 
Thus we can write the entire action as g^ f{<j) where / has 
no g dependence to leading order and the saddle point of 
/ indeed gives exact answers as g — > oo. I will now show 
this just for the quadratic term in the Tr In. It will then 
be clear how it works for higher terms. 

At a schematic level we can write 

Tr\n{D + ga) = Tr{\nD + \n{l + gD^^ a) 

= C^gTrB-^a-^TrD-^GD-^aih.l) 

where C is some constant independent of a. When all the 
indices and details of the Landau interaction are filled for 
this problem and an lo integral is done, we find (dropping 
constants) the following quadratic term 

QT — ^{a\ cos mB cos mO'+crl sin mO sin mO') 



a/3kk' 



-£/3 



.(C(k)<^Mk)</>Mk')0a(k')) 



(5.8) 

where 9 and 9' refer to the angles of k and k' and 
and np are the Fermi factors, and ai and a2 refer to two 
components of the a field. ^ 



^ We need two components because n(k)n(k') cos m(6 — 6') = 
n(k)n(k')(cos cos m6' -|- sinmffsinmS') is a sum of two 
terms and each needs to be factorized with its own Hubbard- 
Stratonovic field. We will refer to the two-component vector as 



We will replace this term by its ensemble average, 
since deviations from the average are down by an ex- 
tra power of 1/g and th us ignorable as g ^ oo 
( Murthv and Mathud . [2003 ) . As for the average, recall 
that to leading order in 1/g, averages of energy levels 
and wavefunctions factorize. Using 



(C(ki)0^(k2)^^(k3)</>Jk4)) 



(^kiki ^kaka 



9 9 

(5.9) 

for our case where ki = k2 = k and ka = k4 = k' and 
summing over the g values of k and k' we find this average 
over wavefunctions goes as l/g. (We are not getting into 
details like a factor of 1/2 coming from the cos^ m9 and 
sin^ m9 in the momentum sums.) 

As for the sum over energy denominators, consider one 
of the two similar possibilities, where = 1 and = 
so that Ea ranges from —gS/2 and while Sp ranges from 
to g5/2. Replacing the sum by an integral over density 
of states 1/5 (valid for the ensemble average) we find an 
integral of the form 



J_ r^^/^ r^^/^ dEdE' _ 5 
7o io E + E' - ~5 

up to constants, so that that 



{QT) 



(5.10) 



(5.11) 



where = a\ + a"^- If we consider higher powers of a 

2 

we again find that they go as . Had we not scaled 

a initially, we would have found that each power of the 
unsealed a was accompanied by a different power of g 
and that no large g limit emerged. 

It can be shown that had we included all Um 's simulta- 
neously, they would not have interfered at the quadratic 
level. Since this is the term that controls symmetry 
breaking, we get the right picture taking just one u,n 
at a time. 

For the experts I mention that since we have a large N 
theory here, it follows as in all large N theories, that the 
one-loop flow and the new flxed point at strong coupling 
are parts of the final theory. However the exact location 
of the critical point cannot be predicted, as pointed out 
to us by Professor Piet Brower. The reason is that the 
Landau couplings defined at a scale El much 

higher than Et (but much smaller than Ep) and their 
fiow till we come down to Et, where our analysis begins, 
is not within the regime we can control. In other words 
we can locate u* in terms of what couplings we begin with 
at Et, but these are the Landau parameters renormalized 
in a nonuniversal way as we come down from Et to Et- 
What is the nature of the state for Um < u^l 
In the strong coupling region <t acquires an expectation 
value in the ground state. The dynamics of the fermions 
is affected by this variable in many ways: quasi-particle 
widths become broad very quickly above the Fermi en- 
ergy, A ( spacings in Vg between successive peaks) has 
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FIG. 8 The phase diagram in the u — - plane. The actual 
phase transition takes place on the line 1/g — 0. However 
it can be perceived at finite g if one is within the l/-shaped 
quantum critical regime. In typical cases, the 1/g axis is 
replaced by the T axis, where T is the temperature. There 
a phase transition of quantum mechanical origin at T = 
influences the quantum critical regime at finite T. 



occasionally very large values and can even be negative, 
and the system behaves li ke one with broken time- 
reversal symmetry if m is oddl Murthv et aH . l2004l ). One 
example of the latter is as follows. Suppose we turn on 
a weak magnetic field B. All quantities- peak positions, 
spacing - will vary linearly in B and not quadratically as 
in a time-reversal invariant system. 

Long ago Pomeranchuk ( Pomeranchuk! . Il958l ) found 
that if a Landau parameter Um of a pure system ex- 
ceeded a certain value, the Fermi surface underwent a 
shape transformation from a circle to a non-rotationally 
invariant form. Recently this transition has rece i ved a 
lot of a t tentio n ( Qganesvan. Kivelson and Fradkinl . l200ll : 
IVarmal Il999fl The transition in question is a disor- 
dered version of the s ame. Details are given in Refs. 
(|Murthv and Shankail . [2003h . (|Murthv et al\ . \2004 . 

Details aside, there is another very interesting point: 
even if the coupling does not take us over to the 
strong-coupling phase, we can see vestiges of the 
critical point and associated critical phenomena. 
This is a general featu r e of r nany qua n tum critical 
poOTfa (|Chakravartv et all . Il988t ISachdevl . Il994l ). i.e., 
points like u^, where as a variable in a Hamiltonian is 
changed, ground state of the system undergoes a phase 
transition (in contrast to transitions wherein tempera- 
ture T is the control parameter). 



How can the cost of adding one particle be negative (after re- 
moving the charging energy)? The answer is that adding a new 
particle sometimes lowers the energy of the collective variable 
which has a life of its own. However, if we turn a blind eye to it 
and attribute all the energy to the single particle excitations, A 
can be negative. 

This result breaks down at exponentially small B, a region we 
can safely ignore since the temperature in any realistic experi- 
ment will be much higer. 



Figure [8] shows what happens in a generic situation. 
On the cc-axis a variable {um in our case ) along which the 
quantum phase transition occurs. Along y is measured 
a new variable, usually temperature T. Let us consider 
that case first. If we move from right to left at some value 
of T, we will first encounter physics of the weak-coupling 
phase determined by the weak-coupling fixed point at the 
origin. Then we cross into the critical fan (delineated by 
the I^-shaped dotted lines), where the physics is con- 
trolled by the quantum critical point. In other words we 
can tell there is a critical point on the x -axis without 
actually traversing it. As we move further to the left, we 
reach the strongly-coupled symmetry-broken phase, with 
a non-zero order parameter. 

In our problem, 1 /g'^ plays the role of T since stands 
in front of the effective action for cr. (Here g also enters 
at a subdominant level inside the action, which makes it 
hard to predict the exact shape of the critical fan.) The 
bottom line is that we can see the critical point at finite 
1/g. In addition one can also raise the actual tempera- 
ture T or bias voltage to see the critical fan. 

Subsequent work has shown, in more familiar exam- 
ples than Landau interactions, that the general pic- 
ture depicted here is true in the large g limit: upon 
adding sufficiently strong interactions the Universal 
Hamiltoni an gives way to other descriptions with broken 
symmetry (|Murth3, \2004) . 

It was mentioned earlier that the critical coupling u* 
(a nonuniversal quantity) cannot be reliably predicted in 
the la rge q limit. It has become c l ear fr om numerical 
work ( Adam. Brouwer and Sharmal . l2003l ) that it coin- 
cides with the bulk coupling for the Pomeranchuk tran- 
sition. In other words, when we cross over to the left u*, 
the size of the order parameter very rapidly grows from 
the mesoscopic scale of order Et to something of order 
the Fermi energy Ep. However the physics in the criti- 
cal fan as well as the weak coupling side is as described 
by the RMT-I-RG analysis. The strong coupling side has 
to be reworked from scratch since the Fermi surface as- 
sumed in the RG that came down to Et is has suffered 
huge deformations (in the scale of Ep). 



VI. SUMMARY AND CONCLUSIONS 

Our goal was to understand the transport properties of 
a quantum dot- an island that electrons could tunnel on 
to and tunnel out of. As one varies the gate voltage Vg be- 
tween the leads and the dot, the conductance G exhibits 
isolated peaks of varying location and height. What we 
wanted was a statistical description of these features as 
exhibited by an ensemble of similar dots. The dot was 
assumed to be so irregular and the classical motion so 
chaotic, that the only conserved quantity was energy. 

To get warmed up, we asked how we would go about 
addressing the problem if electron-electron interactions 
could be neglected. It was seen that peak heights and 
positions could be determined from the wavefunctions 
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and energy levels. These in turn could be determined 
by RMT given just the mean level spacing. In the non- 
interacting theory, P(A), the distribution of differences 
in Vg between successive peaks, was the same as P{S), 
the distribution of spacings between successive levels in 
the dot. 

Actual comparison of these two distributions showed 
a clear disagreement, which represented a failure, not of 
RMT, but of the assumption of non-interacting electrons. 
So interactions had to be included. 

The coulomb interaction implied that adding an elec- 
tron to the dot would cost not just the gap to the next 
empty level, but an additional energy due to repulsion 
by the N electrons already in the dot. This charging 
energy could be accounted for by adding a term uqN'^ 
to the second quantized Hamiltonian. In the presence of 
spin, we also needed to add a term — Jq-S^ to represent 
the exchange interaction. The final result was Hu, the 
universal Hamiltonian. While the success of this model is 
unquestioned at moderate values of interaction strength, 
some arguments for why this had to be the right answer, 
and why other interactions could be neglected because 
their ensemble averages vanished were not persuasive. 
We asked if there was there a better way to understand 
the success of Hu. 

It was pointed out that the KG was such a way since 
it offered an unbiased procedure for determining which 
interactions were really important in deciding low energy 
properties like the ground state and its low energy exci- 
tations. In the RG approach one divided the variables 
into two sets: x, which we cared about, and y, which we 
we did not care about. In our example x was the low 
energy region (near the Fermi energy) and y everything 
else. One then eliminated or integrated out the y's to 
obtain an effective theory of just x, which gave the same 
answers in the x-domain as the original one. In this pro- 
cess some initially very impressive couplings could fade 
into oblivion (irrelevant) while tiny ones could grow in 
size (relevant) and some could remain fixed (marginal). 
In any event we could see which couplings really mat- 
tered. We saw that while the RG concept itself was 
non-perturbative, mode elimination was typically done 
perturbatively in the interaction. 

The application of RG to our problem required a two- 
stage process as developed by Murthy and Mathur. First 
one ignores disorder and finite size and eliminates high 
energy modes outside the Landau band, a region of width 
El measured from the Fermi surface. Here we know from 
past RG work on clean bulk systems that we must in- 
variably end up with the Landau interaction u{9). But 
we are not done yet. Although the Landau scale El is 
much smaller than Ep, the Fermi energy, it does not 
vanish for infinite system size. So Et, which vanishes 
as 1/L, lies even closer to Ep. So we need to renor- 
malize down from El to Et. During this process the 
Landau parameters could renormalize in a way we can- 
not determine. (Though in a clean system Um are strictly 
marginal, once we approach Et and take disorder seri- 



ously, a nonzero flow is guaranteed. ) So one begins at 
Et by writing the (renormalized) Landau interaction in 
the disordered single-particle basis a and eliminating the 
a states to determine the fate of the Landau interactions. 
It was found that only Uq and Jq, the zeroth harmonics 
on the Fermi circle of the Landau interactions (for charge 
and spin densities), survived at the lowest energies if the 
starting value of Um was either positive or not below a 
negative coupling m^, of order unity. Thus the low en- 
ergy fixed point for this range of initial coupling was just 
Hu, the universal Hamiltonian. In addition to providing 
this justification of the emergence of Hu, the calculation 
also suggested that for Um < u*„ the system underwent 
a phase transition in the g oo limit. Before discussing 
the phase transition, let us recall how the calculation of 
the /3-function was done. 

In any interacting field theory with some coupling u, 

one computes a physical quantity like a scattering ampli- 
tude in a power series in u, starting with u itself, followed 
by loop diagrams of increasing complexity. These loops 
involve momentum or energy sums (or integrals) up to 
some cut-off A. Clearly if the sum of all these terms has to 
be independent of A (as is the physical scattering ampli- 
tude) then u itself must become u{A) and vary with A in 
such a way as to keep the series as a whole A-independent. 
Conversely it is possible to determine u{A) by drawing 
diagrams to some order in u making this demand. In our 
problem, the computation of /3(u) = —Adu/dA for any 
one specific dot required knowledge of the wave-functions 
and energy levels lying within Ep. But thanks to self- 
averaging, one could replace the /?- function for the given 
dot by its ensemble average. The zeroes of this /3- function 
are what showed Hu to be a fixed point (at the origin) 
and to be the the critical point for the phase transi- 
tion. 

This clever calculation was nonetheless a weak- 
coupling analysis predicting a phase transition at strong 
coupling. Was the transition real and if so, what was on 
the strong coupling side? It was here that the 1 /N tech- 
nique came in. In the large g limit one could show that 
all physics could be extracted in a saddle point calcula- 
tion of the 1/A^ type for any coupling Um. In contrast 
to theories where there were N equivalent species, here 
we had g fermions with different energies and matrix el- 
ements. However, thanks to disorder self-averaging, one 
could pull out a in front of the action for the Hubbard- 
Stratonovic field cr. The saddle point theory confirmed 
the fixed point nature of Hu for < u„i < oo and the 
transition at u^. Furthermore we could see beyond the 
transition to the other side: here a acquired an aver- 
age (a disordered version of the Pomeranchuk transition 
in clean systems) and produced many attendant conse- 
quences like time-reversal breaking. However the "ex- 
act" critical point Um = u*^ of this calculation was not a 
directly measurable quantity since the saddle point the- 
ory had as its input, not the Landau interaction defined 
at Ep, but what it had evolved into, between El and 
Et- This was a no man's land where disorder was too 
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strong to be ignored, but not strong enough to use RMT 
since we were not within making the flow intractable. 
However clever arguments of Adam et a show that the 
transition occurs at the bulk critical value of Um = 1. 

Since the phase transition occurs only at infinite g (a 
finite system always has finite g and cannot have a tran- 
sition), it might seem that our study of it was academic. 
This is not so, and the reason is the same as in quantum 
critical phenomena. Recall that there, a quantum phase 
transition at T = as a function of some coupling u = u* 
can be felt even at T > inside a F-shaped region called 
the quantum critical region. Here 9^ plays the role of 
1/T as a prefactor in the action. Thus we can see the 
effects of the critical point even at finite g and over a 
wide range of coupling. 

I began with the ominous remark that three obstacles 
are ganged up here: randomness, strong interactions and 
finite size. Yet they ended up being benign: randomness 
and finite size led to a finite Thouless band within which 
we could use RMT and ply our trade, while strong inter- 
actions led to an interesting phase transition. As for the 
KG, it had to be invoked in the first stage as we came 
down to El using the clean system RG to end up with the 
Landau interactions. At this point we could either use 
perturbative self-averaged RG inside Et or better still, 
use the self- averaged 1/N method to solve the model by 
saddle point. 

This colloquium has emphasized what I find most 

beautiful about this problem: the confluence of physi- 
cal complications and the interplay of diverse techniques 
that lead to a solution. Of necessity it has been sparse 
on phenomenology. However, armed with the ideas ex- 
plained here you are ready to remedy this, following any 
number of the excellent references mentioned in the text. 

It will be very interesting if experimentalists unearthed 
the phenomena chronicled here by studying dots with 
strongly interacting electrons, a possibility more readily 
realized than in the bulk since electron density in dots 
can be controlled by gates. Stay tuned for these results. 
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